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We report on new results for the spectrum of quarkonia using a fully relativistic approach on anisotropic lattices 
with quark masses in the range from strange to bottom. A fine temporal discretisation also enables us to resolve 
excitations high above the ground state. In particular we studied the mass dependence and scaling of hybrid 
states. 



1. INTRODUCTION 

The study of heavy quark systems on the lat- 
tice is complicated by the large separation of 
momentum scales which are difficult to accom- 
modate on isotropic lattices. Over the years 
many effective descriptions of QCD have been 
developed for low energies and tested against a 
wealth of experimental data However, the 

(non-)perturbative control of higher dimensional 
operators is very difficult and in practise one 
has to rely on additional approximations. Even 
in non-relativistic bottomonium calculations one 
has observed sizable relativistic corrections to the 
spin structure [^,0 and, more worrying, large 
scaling violations 1| which cannot be controlled 
by taking the continuum limit. We take this 
as our motivation to study heavy quarkonia on 
anisotropic lattices in a fully relativistic frame- 
work. 

The strategy and first results for charmonium 
have already been presented in [0J|]. The ba- 
sic idea is to control large lattice spacing arte- 
facts from the heavy quark mass by adjusting the 
temporal lattice spacing, at, so that m q a t < 1. 
The continuum limit can then be taken at fixed 
anisotropy, £ = a s /a t . The details of our calcula- 
tion are given in Section 2. 

In our study we also employ another advantage 
of anisotropic lattices - a fine temporal resolution 
is ideal to trace the fast exponential fall-off from 
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higher excited states. In particular the gluonic 
excitations have attracted considerable interest as 
they give rise to exotic quantum numbers not al- 
lowed in the quark model. We present our results 
for such exotic hybrid states in Section 3. Section 
4 concludes this report with results from a first 
relativistic bottomonium calculation. 

2. ANISOTROPIC LATTICES 

To study the relativistic propagation of heavy 
quarks it is mandatory to have a fine resolution 
in the temporal lattice direction. To this end we 
employ an anisotropic gluon action: 



s = -P E ^ijW + £ &Pk(x) . (i) 

\x,i>j x,i J 

This is the standard Wilson action written 
in terms of simple plaquettes, P^(x). Here 
(/3,£o) are two bare parameters, which deter- 
mine the gauge coupling and the renormalised 
anisotropy, £ = a s /at, of the quenched lattice. 
The anisotropic gluon action (|l|) is designed to 
be accurate up to 0(a%, af). 

For the heavy quark propagation in the gluon 
background we used the "anisotropic clover" for- 
mulation as first described in The dis- 
cretised form of the continuum Dirac operator, 
Q = m q + tf>, reads 

Q = m + v s Wai + v t Wo7o - 
y [c s (JokFok + c t <JkiF k i\ , 
W„ = V M - (<v/2) 7m A„ . (2) 



This is indeed the most general anisotropic quark 
action including all operators to dimension 5 
(up to field redefinition). The five parameters 
in Eq. || are all related to the quark mass, 
m q , and the gauge coupling as they appear in 
the continuum action. Their classical estimates 
have been given in g. Here we chose m non- 
perturbatively, such that the rest energy of the 
hadron corresponds to its experimental value (e.g. 
M( 3 5") = 3.097 GeV for charmonium). We 
also fix v s — 1 and adjust v t non-perturbatively 
for the mesons to obey a relativistic dispersion 
relation ( c(0) = 1 ): 

E 2 (p) = £ 2 (0)+c 2 (p) p 2 + 0(p 4 )... . (3) 

For the clover coefficients (c s ,Ct) we choose 
their classical estimates and augment this pre- 
scription by tadpole improvement. The tadpole 
coefficient has been determined from the average 
link in Landau gauge: u 0L = 1/3 (f7 M (x)) Land au. 
Any other choice will give the same continuum 
limit, but with this prescription we expect only 
small O(aa) discretisation errors. 

From the quark propagators we construct me- 
son correlators for bound states with quantum 
numbers S(0, 1) x L(0, 1, 2) and, for hybrid states, 
with additional gluonic excitation. Our fun- 
damental bilinears are constructed from two 4- 
spinors, one of 16 T- matrices and additional in- 
sertions of up to 2 lattice derivatives: 

M(x) = q(x)T l A J A k q(x) . (4) 

From those basic operators we can construct 
a vast number of meson states with different 
and definite J PC . For example, the exotic 
quantum numbers 1 ^ can be obtained from 
q(x) eijk Jj eja m A;A m q(x). These simple- 
minded operators can be further improved to give 
states with different projection onto the ground 
state. Here we follow || and employ a com- 
bination of various quark smearings and APE- 
smearing for the gauge links. This allows us to 
extract reliably both the ground state energies 
and their excitations from correlated multi-state 
fits to several channels. We have also checked our 
parameter estimates for consistent fit results from 
different methods and ranges. In order to call a 
fit acceptable we require its Q-value to be bigger 



than 0.1. In Table [l]we list the parameters of our 
simulation. 

A representative compilation of our results for 
the charmonium spectrum is shown in Figure |l| 
for one lattice spacing. Most noticeable is the 
clear level ordering which can be seen for states 
with different orbital/gluon angular momentum. 
To convert our lattice results into dimensionful 
quantities we used the Sommer scale, ro, as de- 
fined in jl0| and calculated very accurately in |ll[] 
for isotropic and anisotropic lattices. It is well- 
known that, without dynamical sea quarks in the 
gluon background, the definition of the lattice 
spacing is ambiguous and one cannot reproduce 
all experimental splittings simultaneously. We 
are however strongly encouraged by the overall 
agreement of this first-principle calculation with 
experimental data (where available). This bodes 
well for the reliability of our quenched predictions 
for states which have yet to be seen by experimen- 
talists. Here we accept the shortcomings of the 
quenched approximation and focus on the scaling 
behaviour instead. 

A detailed study of the spin structure has been 
presented previously (?]||. In the following we 
will analyse the hybrid excitations in more detail. 

3. HYBRID EXCITATIONS 

The exotic hybrid states in Figure [l] are of par- 
ticular interest as they reveal explicitly the gluons 
in the low energy regime of QCD. Here we anal- 
yse the robustness of our results as we change the 
lattice spacing and quark mass in our calculation. 

In Figure ^|we plot our results with fixed renor- 
malised anisotropy, £ = 2, and different spatial 

Table 1 

Simulation Parameters. Using ro to set the scale, 
we adjust the bare quark mass such that the 
1 corresponds to the experimental values in 
ss{cj)),cc(J/y) and bb{T). 
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Figure 1. Quenched charmonium spectrum at 
(/3,£) = (5.7,2). The lattice spacing is set from 
rQ. Horizontal lines show the experimental values. 

lattice spacing, a s . It is apparent that there are 
only small discretisation errors and we estimate 
4.294(71)(200) GeV for the energy of the lowest 
lying exotic charmonium excitation, 1 h , in good 
agreement with other relativistic |T^| and non- 
relativistic || lattice calculations. The second 
error is an estimate for quenching errors which 
presents the biggest uncertainty in our calcula- 
tion. We also expect many non-exotic hybrids in 
this energy region, but they will mix with con- 
ventional states, both experimentally and on the 
lattice. The question whether stable hybrids may 
ultimately be found just below the threshold into 
DiD-dec&y (4.28 GeV) will have to be decided in 
a full dynamical simulation. 

We also find two other exotics, H and 2 H , 
300-500 Mev above the 1 h , as predicted from 
flux tube models jl3]]. Our attempts to measure 
the 0~~ exotic were unsuccessful - presumably 
this excitation is too high to be resolved on our 
lattices @. 

In Figure || we observe a mild mass dependence 
of the exotic spectrum, but the characteristic level 
ordering 1 h < H < 2 H is the same for the 
whole mass range between (f>(ss) and J/^(cc). 

NRQCD simulations have observed a similar 
mass independence of the hybrid levels at even 
higher masses between charmonium and bot- 
tomonium |J. We take this as indication that 
in the quenched theory the hybrid levels are al- 
most completely determined by the gluon dynam- 
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Figure 2. Scaling of exotic charmonium hybrids. 
Expecting mainly 0(a 2 s , a 2 ) errors, we plot our re- 
sults against a 2 (at fixed anisotropy). Within our 
(small) statistical errors we cannot resolve any 
significant discretisation errors for the lowest ly- 
ing exotic 1 '". 

ics (and not the quarks). 

In Figure ^ we also compare our results with 
strange quark mass to a simulation from an 
isotropic lattice |HJ. In contrast to those early 
results we are now able to extract a clear level 
ordering for the lowest lying exotics. We should 
keep in mind that our predictions for the level 
splittings are all within the quenched approxima- 
tion and we should expect deviations of 10 — 20% 
from the real world. 

4. RELATIVISTIC BOTTOMONIUM 

Using the same formalism as described in Sec- 
tion 2, we are also able to perform a relativistic 
calculation of bottomonium where m& ~ 5 GeV. 
To this end we simply adjust the anisotropy to 
have TOf,at < 1. For our initial study we chose 
(/3,£) = (6.1,4) which corresponds to a s ks 0.1 fm 
and at ~ 0.025 fm. The low lying spectrum of 
bottomonium is shown in Figure ||. 

We are very encouraged by those results and 
note the good agreement of the hyperfine split- 
ting (w 30 MeV) with NRQCD simulations at 
this spatial lattice spacing (|^,||. But in con- 
trast to those results in our approach we can con- 
trol scaling violations by taking the continuum 
limit. Work is under way to study the bottomo- 
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Figure 3. Hybrid excitations above the ground 
state. For the lowest lying exotic one can observe 
only a mild mass dependence. 
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Figure 4. Comparison of the high excitations in 
ss with previous work on isotropic lattice Jl5[ . As 
in that study we plot our results normalised with 
respect to the 1 __ . 



nium spectrum on still finer lattices. 

In conclusion, we have demonstrated the ef- 
ficiency of quenched anisotropic lattice calcula- 
tions. A fine temporal resolution is paramount to 
treat heavy quarks in a fully relativistic frame- 
work. Here we also took advantage of a small 
temporal lattice spacing to extract high excita- 
tions in the spectrum for the whole range of quark 
masses from strange to bottom. 
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